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ABSTRACT 

Spectra of Supernovae of type la (SN la) are commonly interpreted as a continuum with absorption 
features. The pseudo equivalent width (pEW) and Doppler shift of absorption features like Sin 
6150A, Can HK 3750A, and the Can PR-triplet 8150A measures which are commonly interpreted as 
the optical thickness and the velocity of a shell of corresponding absorbing material, form the basis 
of common SN la spectral interpretation and classification schemes. 

In this paper, we examine the nature of spectral features in SN la spectra using W7 model spec- 
tra, and show that Sin 6150A and many other features are largely emission dominated instead of 
absorption dominated. We show that apparent absorption features (like Sin 6150A) are frequently 
just coincidental troughs between two (or more) uncorrelated emission features. The pEW measured 
between such emission peaks is then only little related to true strength of the presumed absorption 
feature. It shows how the concepts of "absorption troughs" and "continuum" can be misleading and 
should be used with care. Furthermore, using the same model spectra, we demonstrate for different 
times post explosion how spectral features overlap each other and how together they compose the total 
spectrum. This overlap distorts individual line profiles and affects measured absorption velocities. 

With an explosion model atmosphere that is tuned to match specific observations, the method pre- 
sented in this paper can in principle be used to quantify all these effects and improve the interpretation 
and informative value of observed SN la spectra. 

Subject headings: stars: Supernovae: general - radiative transfer - methods: numerical 
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1. INTRODUCTION 

Supernovae of Type la (SN la) generally sho w very sim- 
ilar spectral and photometric behavior (e.g., Filippcnko 
1997). Because of this interesting homogeneity (and 
the fact that these objects become extremely bright) 
SNe la have proven to be very useful objects to 
study the universe at large distances (|Riess et al.l 119981 : 
iPerlmutter et al]|1999D . Yet despite the general similar- 
ity, there is an even more interesting diversity among 
SNe la, both photometrical and spectral. Whereas pho- 
tometrical diversit y has long been studied (|Phillipslll993t 
IWang et al.l 120031 ) . good spectroscopic data have only 
recently become available in large numbers (CfA SN 
-Prog ram iBlondin et alj |2012t Berkeley SN la Pro gram 
iSilverman et alll2012t) . 

One step towards understanding spectroscopic diver- 
sity is finding similarities and defining classes. A num- 
ber of spectral classification schemes for SN la have been 
proposed in literature. These are based on measure- 
ments performed on conspicuous Si n and Can absorp - 
tion features present in all "normal" ((Branch et allfl993l) 
SN la. These absorption feature measurements are: 1) 
The feature stren gth, exp ressed as pseudo equ ivalent 
width (pEW, see iFolatellil 12001 ([Branch et al.l [20061 ): 
2) The line blue- s hift, expressed as expansion velocity 
(|Wang et all 12009(1 : 3) The rat e of change in blue -shift, 
expressed as velocity gradient (|Benetti et alll200"5l ). 

It is well known that atomic transition lines formed in 
a rapidl y expanding atmosphere, li ke SN la, form P-Cyg 
profiles (jLamers fc Cassinellilll999T ). P-Cyg profiles gen- 
erally show a superposition of an emission wing that is 



more or less symmetric around the transition's rest wave- 
length plus an absorption wing that is blue shifted. But 
the precise shape of P-Cyg profiles in relativistically ex- 
panding atmospheres, and specifically the ratio of emis- 
sion to abso rption, strongly de pends on details of the 
atmosphere (|Hutsemekersl 11993): 1) The radial exten- 
sion of the line forming region; 2) The velocity gradient 
across the line forming region; 3) The wavelength depen- 
dence of the "continuum" and its change across the line 
forming region. 

Parame terized models of SN la atm ospheres, like 
SYNOW (IFishedl^OllBTanch et al.l[2009h or SYNAPPS 
(jThomas et al.l 120111 ). have proven to be very powerful 
and rewarding tools for the interpretation of SN la spec- 
tra. Specifically, these models made it possible to iden- 
tify lines and study the composition of observed SN la. 
However, these models do not account for the geometrical 
extension of the line forming regions nor the velocity gra- 
dients and change of the "continuum" flux across the line 
forming regions involved in shaping the P-Cyg profiles. 
Therefore, such models can not be expected to give accu- 
rate and physically realistic profile shapes, even though 
the large number of free parameters usually enables a 
rather precise reproduction of observed spectra. 

Another fact that complicates the interpretation of SN 
la spectra is the large number of transition lines that 
leave their signature in the spectra in an indistinguish- 
able way (i.e., smeared out by the large exp ansion veloc- 
ities of SN la atmospheres). iKascn (200g) have shown 
that many millions of transitions are needed to calculate 
the SN la opacity accurately. Even though some of these 
lines are stronger than others, the background opacity 
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formed by t he millions o f weaker lines is definitely not 
constant (see[Kasen 2006|) and will leave its signature on 
the profiles of those stronger lines. 

Clearly, the line profiles in SN la spectra contain a large 
amount of information about the detailed atmospheric 
conditions and understanding their nature is the key to 
elevate the use and yields of spectroscopic observations 
to a higher level. In this paper we will demonstrate that 
Sin 6150A and many other features are in fact emission 
dominated and that commonly used measures like the 
pEW and the absorption velocity are little related to the 
actual properties of the these features. 

2. METHODS 

As described in the previous section, a realistic repro- 
duction of P-Cyg profiles requires model atmospheres 
that describe the run of physical conditions over the geo- 
metrical extent of the line forming region. In this paper 
we do not aspire to attain perfect reproductions of ob- 
served P-Cyg profiles but rather focus on understanding 
the general nature of spectral features as they are ob- 
served for SNe la. For this purpose we calculate spectra 
from model atmospheres in which different atomic tran- 
sition lines are removed from the opacity and compare 
them with spectra calculated using the full, unaltered 
opacity. The direct comparison between altered and un- 
altered opacity spectra shows what the signature of the 
removed atomic transition linc(s) on the full spectrum is. 

Because the differences between spectra are often 
small compared to the normal flux level it is impor- 
tant that the model spectra are calculated to very high 
precision. Monte Carlo methods for solving the radia- 
tion transport problem generate spectra in a stati stical 
way, l eading to spectra that are inherently no isy (jLucyl 
Il99a IKasen et all 12005 IKromer fc Simll200l and con- 
sequently are not adequate for this purpose. Instead, we 
use PHOENIX, a st ate-of-the-art, general purpose stellar 
atmos phere code (|Hauschildtl 119921 : lHauschildt fc Baronl 
119991 12004 I2006D . This code solves the special rela- 
tivistic radiation transport problem with very high ac- 
curacy us ing characteristic rays and operator splitting 
methods (Ol son fc KunasJ 119871) . The time evolution 
of SN la atmospheres can be self-consistently calculated 
with PHOE NIX using the radi ation energy balance (REB) 
method (van Rossuml [2012V PHOENIX with the REB 
method can in principle calculate SN la atmospheres in 
NLTE, however the challenging treatment of the huge 
number of lines that play a role for SNe la is work in 
progress (van Rossum, in preparation). Therefore, in 
this paper we use the LTE approximation. 

Determining the atmospheric stratification of SN la 
ejecta from explosion models has been a highly ac- 
tive field of research for decades. Since W7, the 
very successfu l phenomenological ID explosion model by 
iNomoto et al.l ()1984|) , explosion models have continued 
to improve in physical and numerical respe cts (|Khokhlovl 
[19911: lPlewaltaL|[200l: IKasen et al.|[2009l) . Yet all these 
improvements have not necessarily made the model end 
products look more similar to typical SNe la. Therefore, 
over time, W7 has become an established 1-dimensional 
(ID) reference model that is known to reproduce quali- 
tatively the observed spectral properties of SNe la, from 
a few days to months afte r the explosion (IKasen et al.l 
I2006t IKromer fc Siml 120091: Ivan Rossuml 12012ft . making 



this a good starting point for studying the nature of spec- 
tral features. 

Note that with the limitations of W7 (e.g. being ID 
and not specifically tuned to match one certain SN la) it 
can not be expected to perfectly reproduce observations. 
For this reason, we will also not compare res ults to spec- 
tra of individual objects but rather to the Hsi ao et al.l 
(|2007f ) (hereafter Hsiao07) spectral templates. 

3. RESULTS 

We examine the signature of individual lines (or groups 
of lines) by excluding these lines from the opacity and 
comparing the resulting spectrum to the full, unmod- 
ified spectrum. Figure [T] shows the W7 model spec- 
trum compared to the Hsiao07 templates at four dif- 
ferent times post explosion (pe). The peak template 
luminosity is aligned with the peak model luminosity 
at day 18pe. A single scaling factor is applied to the 
unit-less template flux (simultaneously for all times) to 
globally match the model flux. Apart from these two 
tightly constrained parameters, there is no freedom in 
the comparison. The colored patches indicate the differ- 
ence (Delta) between spectra calculated with the lines of 
a certain chemical element removed from the opacity and 
the full spectrum (with unmodified opacities). In order 
to visualize all contributions at once the patches (Deltas) 
are stacked (starting with the low-Z elements), emission 
deltas are stacked below the dashed line and absorption 
deltas above. Patches below (above) the dashed full spec- 
trum curve represent features in emission (absorption). 

Figure [2] zooms in on the Sin 6150A feature in Figure 
[U although shown at different times post explosion to 
complement the information. In Figure[3]the same deltas 
are shown as in Figure Q] but now plotted relative to the 
flux maximum (of the full spectrum) instead of stacked 
on top of the full spectrum and on top of e ach other. As 
a reference, it also shows the lBessell (|1990f ) UBVRI filter 
functions, arbitrarily normalized. 

The effects of line blending on the absorption velocity 
of the Sin 6150A trough are demonstrated in Figure [H 
Figure [5] presents the true equivalent width (EW) val- 
ues for the isolated Sin 6355A, Can IRT, and Can HK 
features, calculated versus the true continuum. 

3.1. Absorption features 

Sin 6150A absorption — On day lOpe, the absorption 
trough around 6150A is dominated by Sin (see also Fig- 
ure[2|). On day 20pe, the Sin feature has mostly changed 
to emission and is significantly supported by Co emission. 
This trend continues on day 30pe, before Fe emission (at 
nearly the same wavelength) completely takes over at day 
40pe. 

Sn "W" absorption — On day lOpe, the absorption 
trough around 5700A is clearly dominated by S II. On 
day 20pe, it is heavily blended with Co and Fe emission, 
and on day 30pe S II has vanished. 

Mgii absorption — On day lOpe, the Mgll absorption 
around 4300 A is blended with Si, S, and Fe absorption. 
On day 20pe, Fe absorption and Co emission dominate 
Mgii, and on day 30pe Mgii has vanished. 



3 



. Ie47 



w7, Phx REB LTE 



u 5 



(d 4 - 




4.0 



le46 



w7, Phx REB LTE 



3000 4000 5000 6000 7000 8000 9000 10000 
Wavelength [Ang] 

Figure 1. W7 model spectra (dashed) are compared to the Hsiao07 
SN la spectral templates (solid) at different times post explosion 
(pe). The comparison uses only one free parameter to simultane- 
ously match the template fluxes. The colored patches represent the 
signature of all atomic transitions of a single chemical element for 
each color. The signatures are individually calculated and subse- 
quently stacked on top of the original spectrum. 

The colored areas above (below) the dashed curve show the absorp- 
tive (emissive) contributions. Apparently, SN la are largely emission 
dominated longwards of 5500Aand absorption features are generally 
much smaller than a pEW evaluation would suggest. 
The famous Sill 6355A "absorption" feature is poorly described by 
the model, especially at day 30pe. The W7 model (and other models, 
see text) seem to miss an important emission feature around 5800A. 
More care should be taken with the interpretation and designation 
of the 6150A trough. 
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Figure 2. Zoom in on the Sill 6150A feature of Figure[T] although 

for different times post explosion: day 6, 8, 12, and 16pe. The 

spectral template flux level is omitted in this plot. 

The true absorption signature of this famous Sill feature is much 

smaller than the total depth of the trough, and decreases over time. 

In these early epochs, the emission wing on the blue side of the 6150A 

trough is mainly affected by O, Si, S, Fe, and Co. The emission wing 

on the red side sees contributions from C, Mg, Si, and Co. 

The wavelength of the trough minimum depends not only on the 

velocity of the Sill 6355 A line but also on how other features are 

blended in (including emission from the Sin 5972A line). 
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Figure 3. Spectral features of different chemical elements nor- 
malized to the unmodified flux peak level in the W7 model (see 
also Figure [TJ . Emission and absorption features are plotted as 
positive an d negative del tas, respectively. The dashed curves show 
the UBVRI Bcssell 1990 band filter functions. 

The famous Sill 6150A and Call HK features, and many others, 
show emission wings that are stronger than their absorption coun- 
terparts. pEW definitions that are based on measuring the depth 
of presumed feature troughs may misrepresent the true feature 
strength. 

The Call HK feature (like many others) does not present itself iso- 
lated but is contaminated with contributions from Si, Co, Ti, and 
S. This complicates both the feature strength and velocity analysis. 
In all four epochs shown here, the B-band (more than other bands) 
contains contributions from many different elements. This possi- 
bly makes the flux in the B-band more robust against signature 
strength (abundance) variations between individual SNe la. 
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Days post explosion 

Figure 4. The absorption velocity of the Sill 6150A feature mea- 
sured in three ways: 1) in the model spectrum calculated with full, 
unmodified opacities (black curve), 2) in the Si curve of the Delta 
plot in Figure|3] (dashed pink curve), and 3) in a Delta plot created 
by removing the Sin 6355 A feature only (solid pink curve). 
The solid pink curve is the most ideal measure of Si absorption 
velocity as it describes an isolated single feature. The dashed pink 
curve shows that the emission wing of the neighboring Si II 5972A 
line (the strength of this line is correlated to the 6355A line) re- 
duces the apparent absorption velocity, although the effect is small. 
After day 28pe, the Si II trough disappears in the model spectra. 
The black curve does not follow the observed, gradually declining 
velocity evolutions because of the 5800A emission feature missing 
in the model. 

The shape of the Si line profile is independent of the presence or 
absence of that feature, so that the solid pink curve still shows a 
realistic run of the true (as opposed to observed) Si absorption 
velocity. 
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Days post explosion 

Figure 5. The evolution of the true equivalent widths (EW) of 
three spectral features over time. Each of the three features has 
a positive and a negative curve corresponding to absorption and 
emission, respectively (note the flipped vertical axis). 
The EWs are computed using the Deltas from individual lines, i.e. 
the line profiles are evaluated versus the true run of the continuum 
flux level over each line profile. 

Both the Sin 6355A and the Call IRT line profiles are emission 
dominated after the first ~10 days pe. The Call HK line profile is 
absorption dominated before day 44pe. 
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Can absorption — The Can HK absorption around 
3750A at day lOpe sees contributions from Si and S ab- 
sorption and Ti emission (best seen in Figure [3]). On day 
20pe and 30pe, the Can HK trough is significantly deep- 
ened by Si absorption, obviously affecting its pEW (note 
that this problem does not affect the true EW values 
plotted in Figure [5]) . This also changes the absorption 
velocity. 

The Can IRT feature appears mainly in emission. 

Fen 5000 A blend — The wide absorption trough around 
5000A is commonly identified as the "Fen blend" (e.g., 
(|Folatellil 12001 IHachinger et al.|[200l V Figure □ and [3] 
show, however, that the Fe is weaker than contributions 
from othere elements at day 10, 20, and 30pe. Only as 
late as day 40pe, this feature is really dominated by Fc 
absorption. 

4. DISCUSSION AND CONCLUSIONS 
4.1. Absence of an important emission peak 

Generally, W7 provides a reasonably good fit to the 
Hsiao07 templates. There is, however, one very seri- 
ous mismatch. The model does not describe the 5800A 
emission peak observed on day 20, 30, and 40pe. (Note 
that this finding is in agreement with results from W7 
calculated with other radiation transp ort codes (e.g., 
IKasen et al.ll200l 1 romcr fc Siml[2009h .^ This is a se- 
rious problem since this emission peak is responsible for 
the blue side of the Sin 6150A trough and thus affects 
the Sin pEW measurements. 

What makes this problem even more serious is that 
other SN la explosion m odels do not descri be this emis- 
sion peak either (e.g., IKasen et al.l 120091 Figure 2b; 
lR5pke et al.ll2012l Figure 3). There are two possible rea- 
sons for this emission peak to be absent from the models: 
1) the chemical element (s) responsible for this observed 
peak is under-represented in the explosion models, 2) 
the opacity data for the elements present in the models 
is missing one or more important transition linetQ. Note 
that NLTE effects do affect the shape of SN la spectra 
(jBaron et al.llTT)!)(i: van Rossum, in preparation) but are 
unlikely to explain discrepancies as big as this missing 
emission peak. 

Since the origin of this observed emission line is un- 
known and it does shape the blue side of the Sin 6150A 
absorption trough, special care must be taken with the 
interpretation of the pEW of this famous trough. 

4.2. Blended features 

Line features are known to strongly blend in SN la 
spectra but still it is common to speak of identified fea- 
tures, named after the predominant contributor. Figure 
Q] and [3] show the degree of blending in the model spec- 
tra, and the way it changes over time, for each feature 
in its own way. Note that the detailed degree of blend- 
ing depends on the atmospheric conditions, which vary 
between SN la objects. The W7 results presented here 
provide a qualitative picture of the blending effects. A SN 
la model atmosphere tuned to fit a certain observation 

1 In the scope of this work, we have verified the opacity data 
used by PHOENIX against recent atomic data from R.L. Kurucz 
(http:/ /kurucz. harvard.edu/ ). 



would allow to quantify the degrees of blending using 
the Delta method presented in this paper. This would be 
a very effective way to increase the informative value of 
observed SN la spectra. 

4.3. The Sin 6 150 A feature 

Since this is the most widely used feature in the char- 
acterization and classification of SN la spectra it deserves 
extra discussion. 

The absorption velocity of the Sin 6355A line versus 
the 6150A trough is demonstrated in Figure |H It com- 
pares the velocity of isolated Si with the value measured 
in the full spectrum. At early times, the Sin 6150A 
trough is filled from the blue side by emission from the 
next bluer Sin line leading to lower velocities, although 
the effect is small (a few percent). 

Figure [5] shows that from day 26pe on the Sin emission 
gradually declines. At the same time the peak is seam- 
lessly taken over by Fe n emission at a similar but slightly 
higher wavelength. This not only affects the pEW but 
also makes the Sin velocity at late times appear lower 
(and its gradient higher) than it really is. 

The run of the full spectrum curve in Figure [5] does not 
match the observed, gradually declining behavior (e.g., 
iBlondin et aHl201~l Figure 15). Also, the 5800A emis- 
sion feature, that is missing in the model (see previous 
subsection), causes the Sin absorption trough to disap- 
pear after day 28pe. This missing feature apparently 
plays a key role in run of the absorption velocity. 

For these rea sons, using the Si l l 6150A velocity at 
late times (e.g., iHowell et all 120061 : IBlondin et aLll2012[ ) 
is potentially problematic. The same pr oblems might 
affect the velocity gradient definition in iBenetti et al.l 
(|2005f h The "improved" velocity gradient definition of 
IBlondin et al.l (|2012t l limits using Sin 6150A velocities 
to day 10 post maximum brightness (approx. day 28pe 
for W7) and thus partly avoids these problems. 

4.4. The nature of spectral features in SN la spectra 

In this subsection we explain how the concepts of "ab- 
sorption troughs" and "continuum" can be misleading 
and should be used with care. 

Misinterpretation of "absorption troughs" — The identified 
features commonly used for interpretation and classifica- 
tion of SN la spectra describe the apparent "absorption 
troughs" , with apparent "continuum" in between of those 
troughs. However, many of the features in SN la spec- 
tra are predominantly in emission. Figure [3] shows that 
for wavelengths larger than 5400A the spectral features 
generally are emission dominated (except possibly Si n 
at very early times), and on day 20pe (close to peak lu- 
minosity, day 18pe) emission contributes significantly to 
the flux level down to wavelengths as low as 3800A. The 
widths and depths of these "absorption troughs" largely 
depend on the strengths and separation of the enclosing 
emission features. 

Misinterpretation of "continuum flux level" — Estimating a 
realistic pseudo continuum from a SN la spectrum is not 
trivial. But since line emission is strong, a "maximum" 
continuum level that connects observed flux peaks is def- 
initely unrealistic for many features at most evolution 
times. 
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Redefinition of SN la feature strength — A SN la model 
atmosphere tuned to fit a certain observation would in 
principle allow to determine the pseudo continuum for 
any spectral feature using the method presented in this 
paper.Yet such atmospheres are not generally available. 

Therefore, instead of the common pEW method, we 
suggest to measure the strength of spectral features us- 
ing the height difference of P-Cyg profiles (at least in 
the region where emission is strong in SN la spectra, see 
above). The stenght is defined as difference in flux levels 
at the trough minimum and the peak maximum on the 
red side of the trough normalized to the mean of these 
levels. The emission peak on the blue side of the trough, 
the strength of which is unrelated to the strength of the 
feature of interest, can thus be excluded, potentially giv- 
ing a more accurate strength measure. Note that this 
new definition is a dimensionless quantity, not taking 
into account the width of the feature, which is loosely 
related to the absorption velocity. Thus the strength 
and the velocity measures become potentially more inde- 
pendent. Application of this redefined strength measure 
to large SN la spectral datasets will show whether this 
measure improves correlations to other spectroscopic or 
photometric properties, and how it will affect classifica- 
tion schemes based on line strengths. 

In future work we will take the method presented in 
this paper to a higher level by using model atmospheres 
that are tuned to match observed SN la spectra. This 
will be parametrized atmospheres, based on W7 and 
more modern explosion models. This way, the details of 
the atmospheric structure can be directly inferred from 
observed SN la spectra in a physically self-consistent way. 
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